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ABSTRACT 

The axion is just one from a general class of new particles - called Light Pseu- 
' doscalar Bosons (LPBs) - predicted by many realistic extensions of the Standard 

t~^. , Model. We offer a somewhat pedagogical review of their main properties, with 

particular emphasis on the effects they induce in a light beam travelling in an 
external magnetic field, like photon-LPB oscillations, birefringence and dichro- 
ism. Moreover, we discuss a new strategy whereby LPBs can be discovered by 
, high-precision observations of certain binary neutron-star systems. Basically, 

in a double pulsar seen almost edge-on, photon-LPB oscillations can give rise 
to a characteristic attenuation pattern of the light beam emitted by one of the 
pulsars when it goes through the magnetosphere of the companion. Depending 
on the actual values of the LPB mass and its two- photon coupling constant, the 
t— > ■ effect can be seen in the 7-ray band with the upcoming GLAST mission. We 

, also shown that this method provides a remarkable cross-check for the recent 

claim by the PVLAS collaboration about the existence of a new LPB. 



1. Introduction 

Light pseudoscalar bosons (LPBs) are predicted by many realistic extensions of 
the Standard Model and have attracted considerable interest in the last few yeari^ 1 . 
Besides than in four-dimensional mo delP, they naturally arise in the context of 
compactified Kaluza-Klein theories^ 1 as well as in superstring theories®. 

Certainly the most well-known example of LPB is the axion, namely the pseudo- 
Goldstone boson associated with the U(l) Peccei-Quinn global symmetry invented to 
solve the strong CP-problen^ . 

As a rule, LPBs are very light and extremely weakly coupled to ordinary matter, 
so that they are quite elusive in accelerator experiments. Still, their characteristic 
two-photon coupling gives rise to photon-LPB conversion effects in the presence of 



a Talk given by M. R. at the "Third International Workshop on Neutrino Oscillations in Venice" 
(to appear in the Proceedings). 



external electromagnetic fields. This phenomenon - quite similar in nature to flavour 
oscillations for massive neutrinos - offers the unique opportunity to detect LPBs by 
means of high-precision non-accelerat or e xperiments. 
In fact, the PVLAS collaboration^ 

has recently reported positive evidence in 
favour of a LPB with mass m ~ 1.0 • 10 _3 eV and two-photon inverse coupling con- 
stanlQ M ~ 3.8 • 10 5 GeV. As we will see, this discovery would definitely signal 
the existence of new physics at low energy which goes far beyond LPBs, and so an 
independent confirmation of the PVLAS claim looks compelling. 

Besides, the two-photon coupling entails profound implications for LPBs in astro- 
physics and cosmologj^ 1 , and - depending on their mass and interaction strength - 
they can also be a natural candidate for nonbaryonic dark matter in the Universe. 

Our aim is twofold. We review the present status of LPB physics, with partic- 
ular attention to the phenomenology arising from their two-photon coupling. After 
this rather general discussion, we address a specific astrophysical implication of the 
photon-LPB conversion mechanism, which leads to a novel strategy whereby LPBs can 
be discovered by high-precision observations of certain binary neutron-star systems®. 
Basically, in a double pulsar seen almost edge-on, photon-LPB oscillations can give 
rise to a characteristic attenuation pattern of the light beam emitted by one of the 
pulsars when it traverses the magnetosphere of the companion. Depending on the 
actual values of the LPB mass and its two-photon coupling constant, the effect can 
be seen in the 7-ray band with the upcoming GLAST mission. We also shown that 
this method provides a remarkable cross-check for the above PVLAS claim. 



2. QED vacuum effects 

In order to better appreciate the physical phenomena through which LPBs can 
be detected, it is convenient to discuss first pure QED effects on the propagation of 
light in a magnetized vacuum. Indeed, we shall see later that qualitatively identical 
effects are brought about by LPBs. Moreover, in order to disentangle the LPB signal 
from the observed data an estimate of the QED background is necessary. 

Our starting point is the lagrangian 
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where the second term on the r.h.s. is the Heisenberg-Euler effective lagrangiarP^ , 
which describes the one-loop corrections to classical electrodynamics - due to the 
exchange of virtual fermions - for photon frequency us much smaller than the electron 
mass m e . Here, a is the fine-structure constant, F^ u is the usual electromagnetic 
field strength, F^ is its dual and natural Lorentz-Heaviside units with h — c — 1 are 
employed throughout. 



3 See eq. (gj) below. 



We assume that an external homogeneous magnetic field B is present, whose 
direction makes an angle 9 with the momentum k of a propagating photon. We 
also suppose that B is much smaller than the critical magnetic field B CT = m 2 e /e ~ 
4.41 • 10 13 G. 

We investigate the effects of B at one-loop level, to lowest nontrivial order in the 
coupling constant. Because of Furry's theorem, diagrams with an odd number of 
photon vertices vanish. So, the triangle diagram with one photon line representing 
the external magnetic field vanishes, and the leading contribution to the photon 
propagator is given by the box diagram with two external-field lines. Clearly this 
contribution is of order a 2 B 2 , but one factor of a can be absorbed into B 2 T so that 
it actually goes like a(B/B cr ) 2 . An explicit calculation shows that the propagation 
eigenstates corresponding to lagrangian Ci are photons with linear polarization either 
parallel or normal to the plane defined by the vectors B and k: they will be denoted 
by 7|| and j±, respectivehyB Owing to the box diagram under consideration, the 
velocity of 711 and j± in vacuo differs from c, and can be conveniently expressed in 
terms of the corresponding indices of refraction. Explicitly, we hav(Pl 
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Thus, we see that the magnetized QED vacuum produces a selective change in the 
velocity of light depending on its polarization, a phenomenon called birefringence in 
analogy with what happens in an anisotropic optical medium. Notice however that 
here birefringence is achromatic, since nn and n± are independent of u. 

As a matter of fact, the above box diagram is relevant also for a different process. 
For, suppose that we replace an external-field line with a propagating-photon line. 
Then what we learn is that in the magnetized QED vacuum the phenomenon of photon 
splitting takes place^D. Even though this argument leads to a correct conclusion, it 
gives a wrong feeling about the orders of magnitudes involved. Indeed, it turns out 
that the diagram in question vanishes in the limit of no vacuum dispersion, namely 
for u — k. Although the previous discussion has shown that vacuum is necessarily 
a dispersive medium in the presence of a magnetic field, an explicit calculation^^ 
reveals that the leading contribution to the photon splitting amplitude does not come 
from the box diagram but rather from the exagon diagram with three external-field 
lines, so that it goes like a 3//2 (B/B cr ) 3 . 

While the significance of birefringence is clear, further discussion is required in or- 
der to bring out the physical relevance of photon splitting. Similarly to what happens 



c We follow the convention to define the photon polarization in terms of the direction of the 
electric field (no general consensus exists on this point). 



for birefringence, the amplitude for photon splitting depends on the polarization of 
the initial photons. A systematic analysis shows that this amplitude vanishes for 711 
photons, so that these photons do not split. Furthermore, it turns out that ■y± pho- 
tons split predominantly into 711 photons^. Hence, we see that the magnetized QED 
vacuum also gives rise to a selective absorption of light depending on its polarization^. 
Thanks to the analogy with what happens in an optical medium, this phenomenon 
is called dichroism. In particular, an unpolarized light beam becomes almost linearly 
polarized in the plane defined by the vectors B and k, owing to photon splitting. 

Actually, the foregoing considerations entail that in the QED vacuum dichroism 
is suppressed by a factor a 1//2 (B/B CT ) relative to birefringence. This fact makes QED 
dichroism totally unobservable in terrestrial laboratories, but it becomes important 
in those astrophysical situations where extreme magnetic fields are known to exist, 
like near the surface of pulsari^. 

It is illuminating to illustrate the above effects in the case of a monochromatic 
light beam which travels in a magnetized vacuum along the z-direction. The beam 
is supposed to be linearly polarized initially (at an angle (p 7^ with respect to the 
plane defined by B and k). 

• Because of birefringence, the two modes ■y± and j\\ propagate with different 
speeds. Therefore at any finite distance from the source, the beam polarization 
turns out to be elliptical. Even more explicitly, as the light beam moves along 
the z-axis, its electric field changes both its direction and its magnitude so as to 
trace a spiral (around the z-axis) with elliptical section. After each 2n rotation, 
a different projected ellipse is singled out in the plane perpendicular to k. Still 
- as long as birefringence alone is concerned - all such ellipses have parallel 
major axis, which is just an elementary manifestation of the composition of two 
harmonic motions along orthogonal directions. 

• Imagine now that only photon splitting is present. Then the resulting dichroism 
depletes the 7_i_ mode, so that the electric field in the 7_i_ mode gets reduced while 
that in the 711 mode gets increased. Geometrically, the net result is a rotation 
of the electric field of the beam. 

As a consequence, when a light beam produced with linear polarization travels in a 
magnetized vacuum an elliptical polarization shows up, with the major axis of the 
ellipses rotated with respect to the initial polarization. 

So far, we have supposed that the external magnetic field is homogeneous. When 
this is not the case, the above discussion is not altered in a qualitative way. Yet, 
a novel effect shows up. Indeed, a varying index of refraction can produce quantum 



d We suppose that u> < 2m e , which ensures that photo-pair production is kinematically forbidden 
(otherwise this process would dominate light absorption). 



vacuum lensing^. This phenomenon is quite similar in nature to optical and gravi- 
tational lensing, apart from the fact that it depends on the light polarization. Observe 
however that QED quantum vacuum lensing is achromatic. 



3. Light Pseudoscalar Bosons 

As we said, LPBs are present in many realistic extensions of the Standard Model 
and their implications for elementary-par ticle ph enomenology, astrophysics and cos- 
mology have been thoroughly investigated^^®. 

3.1. General properties 

A generic feature of LPBs is a CP-conserving two-photon coupling, so that they 
are described by the lagrangian 

A = \ W -\^ 2 ^-^ P 1 " F» v , (4) 

where <fi denotes the LPB field and M is the two-photon inverse coupling constant 
with the dimension of an energjjf). 

A crucial consequence of £2 is that the propagation eigenstates of the photon-LPB 
system differ from the corresponding interaction eigenstates. Hence interconversion 
takes place, much in the same way as it happens for massive neutrinos of different 
flavouri^ 1 . However, since the mixing term in £2 involves "two photons", one of 
them must correspond to an external field. So, photon-LPB transitions occur only in 
the presence of external electromagnetic fields. Observe that such an external field is 
also needed to compensate the photon-LPB spin mismatch. 

As far as a generic LPB is concerned, its mass m and two-photon inverse coupling 
constant M are regarded as independent parameters, which are constrained by the 
known phenomenology (more about this, later). As an orientation, we shall assume 
m < 1 eV in the following. 

The axion naturally fits within this context. Its properties are fully specified in 
terms of the scale f a at which the U(l) Peccei-Quinn global symmetry is sponta- 
neously broken by the vacuum expectation value of a singlet Higgs field. Explicitly, 
the axion mass is 

m .0.6p^X)eV, (5) 
while its two-photon inverse coupling constant is 

M ^ 1 - 2 ' 1010fc (lo4v) GeV ' < 6 > 



e We stress that the following considerations apply also to a scalar light boson, provided the 
F^ v 4> coupling is replaced by F^ v F^ v cp. 



where the constant k is a model-dependent parameter roughly of order one {k — 1 in 
the DFSZ mo del)ff3 Hence, the axion is characterized by the following mass-coupling 
relation 

m ^ .7.*p^)eV. (7) 

Astrophysical and cosmological considerations!^* suggest 10 7 GeV < f a < 10 13 GeV, 
which - thanks to eqs. ([5]) and - entails in turn 10 10 GeV < M < 10 16 GeV and 
10~ 6 eV < m < 1 eV, where we have taken k ~ 1. 

3.2. Photon-LPB oscillations 

At this point, it is appropriate to give a closer look at the phenomenon of photon- 
LPB oscillations. 

Manifestly, they are described by the second- order coupled Klein-Gordon and 
Maxwell equations dictated by lagrangian Ci + C 2 . Quite often - due to the very 
small LPB mass - one is interested in the regime where the photon/LPB energy uj is 
much greater than the LPB mass m itself. In such a situation, the short-wavelength 
approximation can be successfully applied and turns the above wave equation into 
a first-order propagation equation. More specifically, we consider a monochromatic 
light beam travelling along the ^-direction in the presence of an arbitrary magnetic 
field B. Accordingly, the propagation equation takes the following fornA-^ 

i^\1>(z)) = M\1>{z)) , (8) 

with 

|</>(z)> = c x (z) \x) + c y {z) \y) + a,(z) |0) , (9) 

where \x) and \y) are the two photon linear polarization states along the x and y axis, 
respectively, and \<f>) denotes the LPB state. In the {\x), \y), |0)} basis, the mixing 
matrix M. is given by 

B X /2M \ 

M = A yx u + A yy B y /2M . (10) 
V B X /2M B y /2M uj-m 2 /2uj ) 

While the terms appearing in the third row and column of M. have an evident physical 
meaning, the A-terms require some explanation. Generally speaking, they reflect both 
the properties of the medium and the QED vacuum effects addressed in Section 2. 
Off-diagonal A-terms directly mix the photon polarization states and typically give 
rise to Faraday rotation. 

When the external magnetic field is homogeneous, we can choose the ?/-axis along 
the projection of B perpendicular to the z-axis. Correspondingly we have B x = 0, 



By = BsinO, A x = A±, A y = A\\, and the mixing matrix takes the form 

/ u + A ± A R \ 

M=\ A R u + A\\ Bsin9/2M . (11) 
V Bsm6/2M to-m 2 /2u J 

In general, the diagonal A-terms receive three different contributions, and so we write 

A IU = A K D + Afi + Ag . (12) 

The terms Ajj^ D represent the QED vacuum effects. They follow directly from eqs. 
(ED, (ED and read 



In the presence of a plasma, charge screening effects produce an effective photon mass 
given by the plasma frequency otj = Aira n e / m e with n e denoting the electron density, 
and the resulting contribution is 

A^ = -^ = -^. (15) 

"' 2uj m e u 

Furthermore, the term At? ^ describes the Cotton-Mouton effect, namely the birefrin- 
gence of a fluid due to a transverse magnetic field. Finally, the term A# accounts for 
Faraday rotation. 

In the situations to be discussed below, both Ajr^ and Ar turn out to be irrelevant, 
and so they will be discarded from now on. Consequently, the longitudinal component 
of the magnetic field BcosO disappears from M, A±_ decouples while only An mixes 
with 0. As a result, the mixing matrix M. reduces to the two-dimensional form 

( u) + Ay 5sin#/2M \ 
Mo ~ { Bsm9/2M u-m 2 /2u ) ' 

As is well known, such a matrix can be diagonalized by an orthogonal transformation 
with rotation angle 

6 = - arctg ^— — 17 

2 & VAy +m 2 /2uj) K J 

and - in complete analogy with the case of neutrino oscillations^^ 1 - the probability 
the a 7|i photon will be converted into a LPB after travelling a distance d is 
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where the oscillation wavenumber is 
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so that the oscillation length is just L osc = 2n/ A osc . Notice that -P(7y — > 4>) becomes 
energy-independent when A osc is dominated by the photon-LPB mixing term. 

3. 3. Photon propagation 

A further consequence of the photon-LPB mixing concerns the propagation of 
light. As first pointed out by Maiani, Petronzio and Zavattinf^ 1 , this actually comes 
about in two distinct ways. 

We have seen in Section 2 that the exchange of virtual fermions yields a nontrivial 
contribution to the photon propagator in the presence of an external magnetic field, 
which affects the velocity of light and gives rise to vacuum birefringence. A further 
contribution of the same kind arises from the exchange of virtual LPBs, namely from 
the diagram in which two 770 vertices are joined by the line and two 7 lines represent 
the external field. Again, the corresponding amplitude depends on the polarization 
of the initial photon, and so it is an additional source of vacuum birefringence. As 
before, the photon propagation eigenstates are the linear polarization modes 711 and 
7_i_, but a small complication arises here, due to the fact that the mass eigenstates 
differ from the interaction eigenstates. More specifically, 711 and 7j_ now denote the 
polarization modes of the photon states which diagonalize the mixing matrix A4. We 
suppose for simplicity that the magnetic field is homogeneous. Then we know that 
the 7_i_ mode decouples from 0, and hence its refractive index is still given by eq. (J3]). 
On the other hand, a straightforward calculation^ 1 shows that eq. (j2J) gets presently 
replaced by 



Since we expect a similar result to hold true also when the magnetic field is inhomo- 
geneous, we see that the resulting quantum vacuum lensing becomes chromatic in the 
presence of LPBs. 

Next, we address the effect arising from the production of real LPBs, occurring 
via the 770 vertex with a photon line representing the external field. Because only 
711 photons mix with LPBs, photon-LPB conversion is an additional source of vacuum 
dichroism. 




(20) 



Consider now a monochromatic light beam which travels in a magnetized vacuum 
along the z-direction, and suppose that it is linearly polarized initially (at an angle 
ip 7^ with respect to the plane defined by B and k) . How is its propagation affected 
by the existence of LPBs? We know that birefringence gives rise to an elliptical 
polarization, whereas dichroism produces a rotation of the ellipse's major with respect 
to the initial polarization. Manifestly, these effects are qualitatively identical to those 
arising from the QED magnetized vacuum, but the point to be stressed here is that a 
measurement of both the beam ellipticity and the rotation angle - induced by LPBs 
- permits a complete determination of their two independent parameters, namely m 
and M. Such a circumstance turns out to be of immediate relevance for the discovery 
of LPBs, as we will see in Section 4. 

3.4- Astrophysical constraints 

A crucial implication of the two-photon coupling of LPBs concerns stellar evolu- 
tion. Thermal photons in the inner regions of stars are easily converted into LPBs in 
the fluctuating electromagnetic fields of the stellar plasma. So long as LPBs are fully 
described by lagrangian C 2 , they traverse unimpeded the stellar material and escape 
into free space. Hence, the stellar energy-loss rate gets increased in the presence of 
LPBs. Because of the negative specific heat resulting from virial equilibrium, the 
central temperature also increases, thereby changing the observed properties of stars. 
Yet, current models of stellar evolution are in fairly good agreement with observa- 
tions. Therefore, the two-photon inverse coupling constant M has to be large enough 
to provide a sufficient suppression of unwanted LPB effects. 

This argument has been applied in a quantitative fashion to stars of different 
kindU, with the result 1 ' that typically M > 10 10 GeV. 

Remarkably enough, basically the same conclusion is reached by the CAST ex- 
periment, designed to detect LPBs (axions) emitted by the Sun^ 1 (more about this, 
later) . 

4. Laboratory experiments 

Over the last two decades, various proposals have been put forward to detect ax- 
ions through their two-photon coupling in non-accelerator situations. Generic LPBs 
can also be discovered by these techniques, provided of course that their mass and 
inverse two-photon coupling constant fall into suitable ranges dictated by the exper- 
imental setup. 

A strategy addresses the possibility of converting actually existing axions into 
photons. As suggested by Sikivie^J, Galactic dark matter axions are expected to 
excite a proper TM mode of a tunable microway cavity - permeated by a strong mag- 



We do not need to be more precise about this point here. 



netic field - when resonance takes place, namely when the characteristic frequency of 
a cavity TM proper mode happens to coincide with the axion mass. So far, no posi- 
tive signal has been reported by experiments of this kincPU 1 . Although this method 
can also be applied to axions coming from the Sun, they can be more easily detected 
by back conversion into X-ray photons in a "magnetic telescope"^!. Recently, the 
latter technique has been implemented by the above-mentioned CAST experiment!^ 1 , 
leading to the lower bound M > 0.86 ■ 10 10 GeV for m < 0.02 eV. 

An alternative strategy relies upon virtual LPB effects produced on photon prop- 
agation. As we have seen, an initially /mear/?/-polarized light beam travelling in a 
strong magnetic field is expected to become elliptically polarized, with the ellipse's 
major axis rotated with respect to the initial polarization. Although an induced ellip- 
ticity is also due to the magnetized QED vacuu nfl, Maiani, Petronzio and Zavattini 
have shown that the LPB contribution dominates for the physically interesting region 
of the m — M parameter plane in the case of the axiorA^ 1 . Recalling the conclusion 
reached at the end of Subsection 3.3, an independent determination of m and M 
can be achieved by measuring the axion-induced beam ellipticity and rotation angle. 
Needless to say, also LPBs can be searched for in this way At variance with the pre- 
vious technique, here it is clearly irrelevant whether LPBs are dark matter particles 
or whether they are emitted by stars. 

Using an experimental setup based on the latter method, a few years ago the 
lower bound M > 2.8 • 10 6 GeV was established for m < 10~ 3 eV^^ 1 . Quite recently 
- by exploiting a similar technique - the PVLAS collaboration has reported positive 
evidence for an anomalously large value of the rotation angle in an initially linearly- 
polarized laser beam undergoing multiple reflection in a 5T mag netic fielcPS. In 
addition, also the beam ellipticity has been determined. Assuming that the effect is 
indeed brought about by a LPB, the corresponding physical parameters turn out to 
be in the range 1.0 • 10" 3 eV < m < 1.5 • 10~ 3 eV and 2 • 10 5 GeV < M < 6 • 10 5 GeV. 
In the following, we will adopt the values^' m ~ 1.0 • 10 -3 eV and M ~ 3.8 • 10 5 GeV. 

Manifestly, a look back at eq. (J7J) shows that the LPB in question cannot be 
the axion. Moreover, the quoted value of M violates both the astrophysical bound 
discussed in Subsection 3.4 and the CAST result by about five orders of magnitudes. 
The only way out of this difficulty is evidently to suppose that such LPBs are not 
emitted by stars, which amounts to say that either they are not produced at all or else 
they are confined inside the star interiors in a manner consistent with the observed 
properties. In either case, new physics at energy as low as a few KeV is required. 
This is an exciting possibility which is becoming to be addressed^ . At any rate, the 
need for independent tests of the PVLAS result looks imperative. 

Concerning the latter point, a few options have been contemplated. In the first 
place, an experiment similar to PVLAS can be performed, however with a different 



s No observable rotation is expected to arise in the QED vacuum because dichroism is suppressed 
with respect to birefringence, as pointed out in Section 2. 



magnet and a laser beam with a different frequency, like in the BMV project!^ 1 . A 
somewhat different method exploits the idea of photon regeneration^ in connection 
with a synchrotron X-ray beam produced in a free-electron laser, like the LCLS at 
SLAC and the XFEL at DES^23. Alternatively, the value of M found by PVLAS can 

IOOI 

be tested by data which will become available in the near future at e + e~ colliders^, 
since in these machines the two-photon coupling of a LPB gives rise to events with a 
single photon and missing energy in the final state. 



5. Double-pulsar observations 

A remarkable astrophysical implication of the photon-LPB conversion mechanism 
concerns a binary neutron-star syste nP 

and can be illustrated as follows. 

Consider a binary pulsar seen almost edge on and focus the attention on the 
light beam coming from one of the two stars, to be referred to as "pulsar A" . Sup- 
pose further that both components are neutron stars, so that basically no accretion 
phenomenon perturbs the system. Denote by p the impact parameter of the beam, 
namely its projected distance from the companion, that is to say from "pulsar B" . 
When p is large, nothing interesting is expected to happen. But when p is small, 
the beam goes through the magnetosphere of pulsar B, so that the observed pho- 
tons traverse a region where a strong magnetic field is present. As a consequence, 
some of these photons are expected to be converted into LPBs, thereby producing a 
characteristic attenuation pattern of the beam intensity. 

Realistically, we consider the recently-discovered binary pulsar system J0737- 
3039^', which is particularly well suited for our purposes. Indeed, its orbital plane 
has currently an inclination angle i ~ 87°, and the orbital precession of the system 
allows us to predict that the value % ~ 90° could be attained before 2020. The or- 
bital period is 2 hours and 45 minutes, while the spinning period of pulsar B is 2.77 
seconds. Although the current minimum value of the impact parameter of the light 
beam from pulsar A is somewhat uncertain, an estimate^ 1 yields p ~ 4 ■ 10 3 km and 
this is the value used in our analysis. Depending on the actual values of m and M, 
the beam attenuation can be substantial. We stress that the effect in question does 
not depend on any new physics responsible for the confinement of the LPBs in the 
star interiors. 

Rather than carrying out a general analysis for arbitrary values of m and M, we 
prefer to focus here the attention on the specific values corresponding to the PVLAS 
result, namely m ~ 1.0 • 10~ 3 eV and M ~ 3.8 • 10 5 GeV. 

The geometry of the considered double pulsar is depicted in Fig. 1. We introduce 

h Suppose that a light beam is shone across a magnetic field, so that some photons are converted 
into LPB. If a screen is put on the beam path, photons are completely absorbed, but LPBs are 
not. So, they emerge on the other side of the screen and undergo photon back conversion if another 
magnetic field is present. Therefore, the presence of photo ns on the latter side of the screen would 
be an unambiguous signal of photon-LPB conversiorPSESl 



a (x,y,z) coordinate system fixed with pulsar B. The light beam from pulsar A 
propagates along the z direction and passes pulsar B at impact parameter p. For 
simplicity, we suppose that the rotation axis of pulsar B lies along the y direction 
and that its magnetic dipole moment precedes at 45° about the y axis. Therefore, the 
light beam from pulsar A experiences a time-dependent and inhomogeneous dipolar 
magnetic field B, with typical strength B ~ 10 12 G on the surface of pulsar B. In such 
a situation, the mixing matrix has the general form (JTUl) . Since we are not interested 
in the polarization of the photon beam, we can safely set 

A xy ~ A yx ~ . (21) 

For the same reason, we can take 




on account of eqs. f[T3"l) and (|14p . Plasma effects arise from the presence of free 
charges in the magnetosphere of pulsar B and their contribution is given by eq. (TT~5T) 
with the electron density supplied by the Goldreich- Julian relatiorPH 1 

n e ~ 7 • 10- 2 (|) (~) cm- 3 , (23) 
where P denotes the spinning period of pulsar B. 




Figure 1: Geometry of the model binary pulsar system 

We proceed to solve the resulting eq. (|HJ) numerically in the equivalent exponential 
form 

\^{z)) = Texp (-i £dz'M(z'f) |V(zo)) , (24) 

where T is Dyson's ordering operator. Since the A4 matrix is explicitly z-dependent, 
the propagation is computed iteratively, i.e. in steps of small intervals 5z. For small 



enough 8z, the M. matrix can be considered as constant over a single distance step, 
and the short-propagation iterations can be calculated as 

\ip(z + 5z)} ~ exp(-iM{z)6z)\ip{z)} . (25) 

The propagation is performed from Zq to z mSuX , with the initial condition corresponding 
to a pure photon state. Here, z is fixed by the trajectory of pulsar A and z max is 
chosen in such a way that B x (z)/2M ~ B y (z)/2M ~ for z > z max . The transition 
probability is then 

P{li -0) = |<</#0w))| 2 = |c (^ max )| 2 , (26) 

where j = x (y) refers to the initial photon polarization along the x-axis (y-axis). 
Clearly, the transition probability for an unpolarized photon beam is 

P (7 ^ ) = S2=^ + £C*^*>. (27) 

We plot P(7 — > <p) versus photon energy uj in Fig. 2, for the beam impact 
parameter p = 4 • 10 3 km. 




Figure 2: Transition probability versus photon energy for a trajectory of the light beam with an 
impact parameter p = 4 • 10 3 km. 



We see that photon-LPB conversion turns out to be important for u > 10 MeV, 
namely in the 7-ray band. This is a remarkable conclusion, since pulsar A in J0737- 
3039 is expected to be a 7-ray source. In addition, 7-ray photons propagate totally 
unimpeded in the magnetosphere of pulsar B, and so we do not have to bother 
about further potential sources of beam attenuate Finally, the short-waveleagth 
approximation - upon which eq. (jHJ) is based - now appears to be fully justified. 



'We stress that in the present situation photo-pair production is totally negligible. 



Actually, the behaviour of P(7 — > 0) in Fig. 2 can be understood in an intuitive 
fashion. For, imagine a drastically simplified situation in which the magnetic field 
experienced by the light beam is stationary and homogeneous. As already pointed 
out, a typical value of the magnetic field strength on the surface of a neutron star 
- namely at a distance r ~ 10 km from the centre - is B ~ 10 12 G. Since B is 
dipolar, it goes like ~ r -3 , and so the magnetic field strength at p ~ 4 • 10 3 km 
1.6 ■ 10 4 G. Correspondingly, from eq. ( 1221) we have A^f D ~ ^yy D — 
~ 18 (cj/MeV) eV. Owing to eqs. (fPol) and ff23l) . it is straightforward to see 



m b 

6.7 ■ 10" 



that the term A£r ~ A^f is much smaller than AS,^ D 



yy 



Agf D for 5 



1.6- 10 4 G, 

and so it can be neglected. Moreover, we find B X /2M ~ B y /2M ~ 0.4 • 10~ 12 eV 
and m 2 /2uj ~ 0.5 • 10~ 12 (MeV/u;) eV. Because the terms A XJ/ , A y:r presently get only 
a QED contribution, they should not exceed A^f D ~ A^ D . As a consequence, 
condition (I2T!) appears to be justified. Now, a look back at eqs. (TT7T) and (|T8l) shows 
that mixing effects in A^o should become physically important for B ~ 1. Therefore, 
substantial photon-LPB conversion is expected to take place for u > 1 MeV, which 
explains the threshold behaviour of P(7 — > 0). Furthermore, we know that this 
quantity becomes energy-independent when the photon-LPB mixing term is much 
larger than the others. This situation presently occurs for u ^> 1 MeV, thereby 
explaining also the flat behaviour of P(7 — ■> <p). 

Assuming that the light beam emitted by pulsar A is not polarized, the temporal 
behaviour of the considered effect is best expressed in terms of the total transmission 
A = 1 — P(7 — > 0) of the beam after propagation in the magnetosphere of pulsar B. 
We plot A versus time in Fig. 3, as pulsar A moves in its nearly edge-on orbit. 
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Figure 3: Total transmisson of the gamma photon beam emitted by pulsar A versus time. Inset 
shows the modulation mainly due to the rotation of the magnetic dipole moment of pulsar B. 



Our numerical simulation predicts a strong attenuation of the photon beam up 
j\Vith a magnetic field of this strength, photon splitting is completely irrelevant. 



to 50 % with a time duration of about 200 s. As it is clear from Fig. 3, this effect 
has three different temporal structures. The broad minimum - from —200 s to +200 
s - evidently corresponds to the transit of pulsar A behind pulsar B. The tens-of- 
seconds, symmetric peaks are due to photon-LPB oscillations, depending on the actual 
path through the interaction region with pulsar B. Finally, the highest frequency 
modulation - shown in inset - is due to the rotation of the magnetic dipole moment 
of pulsar B. 

We note that pulsar A also generates a similar effect on the light beam emitted 
by pulsar B. Since the spin period of pulsar A is smaller than the one of pulsar B 
(23 ms), the resulting modulation period will be also smaller. 

Let us finally demonstrate that the effect under consideration is observable with 
the upcoming GLAST mission, which will have a pulsar 7-ray detection sensitivity 
up to 2 orders of magnitude better than any previous telescope^^ 1 . To this end, we 
show in Fig. 4 the region of the m — M parameter plane excluded by the no detection 
of an attenuation A at the 10 % level. 
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Figure 4: Exclusion region in the case that the existence of the attenuation is excluded at 10 % level. 

Such an attenuation is achieved by a 100 photon count during the total integration 
time. For a two weeks observation time, this corresponds to a flux from pulsar A 
of about 2 • 10~ 7 photons / cm 2 /s , which is a reasonable flux according to previous 
observations of several pulsars^ 1 . As a matter of fact, the GLAST sensitivity curves 
allow for a much weaker minimum detectable fluxP^ 1 , down to 1-10~ 10 photons/cm 2 /s. 

Thus, we conclude that the photon-LPB conversion mechanism in the double 
pulsar J0737-3039 really provides a cross-check for the recent PVLAS claim about the 
existence of a new LPB. We stress that our result would remain practically unchanged 



even if we were considering a light scalar 
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